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Three kinds of integrable Kondo problems in one-dimensional extended Hubbard models are 
studied by means of the boundary graded quantum inverse scattering method. The boundary K 
matrices depending on the local moments of the impurities are presented as a nontrivial realization 
of the graded reflection equation algebras acting in a (2sc + l)-dimensional impurity Hilbert space. 
Further, these models are solved using the algebraic Bethe ansatz method and the Bethe ansatz 
equations are obtained. 
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I. INTRODUCTION 



The study of integrable models of correlated electrons with open boundary conditions has been the subject of 
considerable attention 1^-^. Recently it has become apparent that for models on open chains it is possible to obtain 
integrable impurity boundary conditions as operators which need not be expressed in terms of the (super)symmetry 
of the bulk model. A very important application of this procedure is in the context of Kondo; i.e. spin impurities 
in models of correlated electrons. For the case of the supersymmetric t — J model boundary spin-i impurities were 
introduced in | [lO| ] and the resulting model solved by means of the co-ordinate Bethe ansatz method. 

A reformulation of this model in the context of the Quantum Inverse Scattering Method (QISM) was given in |Tl| ] 
demonstrating that the model could be obtained via a family of commuting transfer matrices and thus establishing 
integrability. Central to this approach is the representations of the reflection equation algebras originally introduced by 
Sklyanin . Such a solution guarantees that boundary terms may be applied to any model whose bulk integrability 
is associated with a solution of the Yang-Baxter equation. An interesting observation made in pd] | was that the 
necessary solution of the reflection equation was not regular in the sense that it is not obtained by "dressing" ; i.e. it 
can not be factorized into a product of local monodromy matrices and a c-number matrix. 

By utilizing the underlying algebraic structure it was subsequently shown in that a more general classes of 
integrable t — J models with Kondo impurities exist. These were derived from both gl{2\l) and (7^(3) invariant solutions 
of the Yang-Baxter equation and the solution of the reflection equation was extended to accomodate arbitrary spin s 
impurities situated on the boundaries. Again, the new solutions of the reflection equation are not regular. Moreover 
it was also demonstrated in |13| that the algebraic Bethe ansatz is applicable for these models and explicit solutions 
were given. 

Recently, the work of Frahm and Slavnov has provided a representation theoretic explanation for the existence 
of these non-regular solutions of the relection equation. In essence, such solutions are obtained by suitable projection 
onto a subspace of the impurity Hilbert space for a regular solution. A consequence of this projection method is that 
the remaining (super)symmetry in the new boundary operator on the impurity site corresponds to a subalgebra oof 
the (super)symmetry of the original regular solution. As examples, this was illustrated in JT^ for the case of gl{m) 
impurites coupled to an open gl{n) invaraint chain for m < n and a reproduction of the integrable t — J model with 
Kondo impurities given in [ p^ . 

It is immediately evident in view of these results that integrable spin impurities, being characterized by the simplest 
Lie algebra su{2), can be readily obtained from regular solutions coming from the larger (super)symmetry associated 
with the model in the bulk. In particular, it is possible to obtain integrable boundary Kondo impurity models 
associated with the Lie algebra gZ(4) and superalgebras gl{3\l) and gl{2\2) which we investigate here. In each case, 
the bulk Hamiltonian can be expressed in the form of an extended Hubbard model and thus is worthy of investigation 
in terms of the physical properties that are exhibited. The bulk Hamiltonian associated with the gl{2\2) solution is 
well known from previous works of Essler et. al. p5| ]. However the other two cases give rise to bulk Hamiltonians 
which are apparently new. 

In the next section we introduce the three forms of extended Hubbard models with integrable boundary Kondo 
impurities. Following this we undertake an algebraic Bethe ansatz approach to solve each case. In the last section we 
conclude with some final remarks. 
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II. INTEGRABLE NON-C-NUMBER BOUNDARY /f-MATRICES AND KONDO IMPURITIES IN 
ONE-DIMENSIONAL EXTENDED HUBBARD MODELS 



Let Cj.o- and cj ^ denote fermionic creation and annihilation operators for spin a at site which satisfy the anti- 
commutation relations {c\^,Cj^r} — SijSar, where i,j — 1,2,---,L and tr, r =t, j. We consider the following 
Hamiltonian which describes two impurities coupled to the supersymmetric extended Hubbard open chain of Essler 
et. al. 

L-l 



H = - + H.c.)(l - rij.-<T - nj+i^-a) 

L-l L-l ^ 



i=i j=i 

+ JaSi ■ Sa + Vafli + UaTli-^nii + JfeS^ • Sfc + VfoTT-L + UbUL^riLi, (1) 

where Jq, Va and = a, 6) are the Kondo couphng constants ,the impurity scalar potentials and the boundary 

Hubbard- like interaction constants, respectively; S is the vector spin operator for the conduction electrons; Saia = a, b) 
are the local moments with spin-i located at the left and right ends of the system respectively; rija- is the number 
density operator rijo- = Cj^cja, nj — riji + Uji. 

The supersymmetry algebra underlying the bulk Hamiltonian of this model is gl{2\2). It is quite interesting to note 
that although the introduction of the impurities spoils the supersymmetry, there still remains u{2) ® u(2) symmetry in 
the Hamiltonian (|^) whose representation contains the spin and 77-pairing realizations. As a result, one may add some 

terms like f/X]^=i '^jT'^iiiMS^=i ^X]j=i("iT ~ "■ji) ^'^ ^^'^ Hamiltonian (0), without spoiling the integrability. 

Below we will establish the quantum integrability of the Hamiltonian (|l|) for a special choice of the model parameters 
Jq, T4,,and Ua 

J 2 ^ -t- 2CaSa jj 2Sq, Cq(25q 1) 

+ 1)' " C„( -fl)' Ca{Ca+2Sa + l) 

This is achieved by showing that it can be derived from the (graded) boundary quantum inverse scattering method 
|^,|j.Here we emphasize that a special case of this model, corresponding to Sq = ^, has been studied in p6| |. 
The second choice of couplings which leads to an integrable model is given by 

L-l 

H = - (c] „Cj+i,^ + H.c.)(l - nj- _^ - nj+i,_^) 

L-l L-l g 

L-l L-l 

-2 51 "jM'^j:T('^i+ia"j+i,T - nj+i) - 2 ^ '^j+i,i'^j+i,TK-a"j,T " "j) 

+ JaSi • Sa + VaUi + C/a"-lT"U + -h^L ' Sf, + VbTlL + UbUL^nLl, (3) 

In this case we can introduce integrable Kondo impurities on the boundary by choosing 

8 4cl + 4cc,-Asa{sa + l)-3 _ Acl + 8Cc, - 4Sa{Sa + I) - 5 



(2Cc. + 2Sc + l)(2Ca - 2s„ - 1) (2Ca 2Sa + l)(2Ca - 2s„ - 1) (2Ca + 2s„ + l)(2Ca - 2sa - 1) 

(4) 

A third choice of couplings which leads to an integrable model is 

L-l 

H = - {c]^„Cj+i.a + H.c.)(l - - %+1,-a) 

L-l L-l ^ L-l 

i=i j=i j=i 

-l-JaSi • Sa -f V"a«l + C^a'T-lT"U + ^^fc^i ' ^fa + VfeTT-L + UbTlL^nLi, (5) 
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where integrable Kondo impurities on the boundary are obtained by the choice 

(24 -l)2_4s„(s„ + l) 



J cy. — 



(2Cq + 2Sa + l)(2Ca - 2sc - 1 



(2c„ + 2sa + l)(2c„-2s„-l) 



(2ca + 2sa + l)(2ca - 2s„ - 1) ■ 

(6) 



Let us recall that the Hamiltonian of the ID supersymmetric extended Hubbard model with periodic boundary 
conditions commutes with the transfer matrix, which is the supertrace of the monodromy matrix T(u) 



T{u) = i?oL («)■■• Roi {u) . 

Here the quantum R-matrix R[u) comes from the fundamental representation of gl{2\2) and takes the form 
/ 



(7) 



R{u) 



- 2 
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(8) 



V 

It should be noted that the supertrace is carried out for the auxiliary superspace V . The elements of the supermatrix 
T{u) are the generators of an associative superalgebra A defined by the relations 



i?12(Ml - U2) T (Ui) T iu2) =T (M2) T (Mi)i?i2(wi - U2), 



(9) 



1 2 

where X= X X= I® X for any supermatrix X £ End{V). For later use, we list some useful properties enjoyed 
by the R-matrix: (i) Unitarity: Ri2{u)R2\{~u) — p{u) and (ii) Crossing-unitarity: i?^2^(— U+ l)i?2i^(u) = p{u) with 
p{u),p{u) being some scalar functions. 

In order to describe integrable models on open chains, we introduce two associative superalgebras 71 and defined 
by the R-matrix R{ui — U2) and the relations 



1 22 1 

i?12(wi - U2) T- (Ml)i?2l(ui + U2) T- (U2) =T- (M2)i?12(wi + "2) T- (Ml)i?2l(ui " U2) , (10) 



+U2) rf^ {ui){ [Rf,'{u,+U2) r'^ tI^'- {u2) 

=Tf- {U2){ [ (^1 + «2) {ur)Rxr''{-Ux + U2), (11) 

respectively. Here the supertransposition sta (a = 1, 2) is only carried out in the a-th factor superspace of y ® 1^, 
whereas ist^ denotes the inverse operation of sta- By modifying Sklyanin's arguments |T2[ , one may show that the 
quantities r(w) given by r(w) = str{T+{u)T-{u)) constitute a commutative family, i.e., [t{ui),t{u2)] — 0. 
One can obtain a class of realizations of the superalgebras 7^ and 71 by choosing T± (u) to be the form 

r_(u) = T.{u)f.iu)TZ\-u), Tfiu) = Tl\u)ffiu) {T-\~u)f (12) 

with 
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T„(w) = RoAiiu) ■ ■ ■ i?oi(u), T+{u) = i?oL(u) ■ ■ • i?o,Af+i(w), T±{u) = K±{u), 



(13) 



where K±{u), called boundary K-matrices, are representations of 7± in some representation superspace. 

We now solve ( |l0| ) and for K^{u) and _ft'+(w). For the quantum _R-matrix (^), one may check that the matrix 
K- (u) given by 



/ 1 \ 

10 

A^{u) B^{u) 

V c_(w) y 



(14) 



where 



B-{u) 
C-{u) 
D-(u) 



+ 2U - Acl ~ ^Ca{2Sa + 1) + 4uSg 
{u-2Ca){u-2Ca-'^Sa-2) 

4tiS^ 

(w-2ca)(u-2ca-4s„-2)' 

4uSj 

(u-2c,)(u-2ca-4s„-2)' 

+ 2» ~ Acl ~ 4c, (2.9, + 1) - 4uSg 

[u-2Ca){u-2Ca-ASa-2) 



(15) 



satisfies (|^. Here S± = ± iS^. The matrix X+(w) can be obtained from the isomorphism of the superalgebras 
71 and 7I|_. Indeed, given a solution T_ of (p^, then 7+(u) defined by 



Tf(«) = r_(-u) 



(16) 



is a solution of ( [Til ) . The proof follows from some algebraic computations upon substituting (|T^ into ( [Til ) and making 
use of the properties of the R-matrix . Therefore, one may choose the boundary matrix _ftr_|_(u) as 



K+{u) 



/ 1 \ 
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A+{u) B+{u) 

\ C+{u) D+iu) j 



(17) 



with 



A^{u) 
B+{u) 
C+{u) 



- 2u - Acf - Acb{2sb - 1) + : 



{u - 2cb + 2){u- 2cb - Asb) 

4wSr 



{u - 2cb + 2){u- 2cb - 4s6) 
4^S+ 

(u-2cb + 2)(u-2c6-4sfc)' 
1*2 - 2m - Acl - 4c6(2sf, - 1) + 8sb - 4uS^ 
(u - 2cf, + 2)(u - 2cb - Asb) 



(18) 



Now it can be shown that Hamiltonian (|i|) is related to the second derivative of the boundary transfer matrix t(u) 
with respect to the spectral parameter u at m = (up to an unimportant additive constant) 



H 



r"(0) 



L-l 



A{V + 2W) 



2{V + 2W) 



stro K+ (O)Glo 



+2 stro ( A''+ (O)iJ'^ 



LO 



Stro ( /^+ (0) (i?fo)' 



(19) 



with 



1 d 

2du 



PR{u) 



A 



where P denotes the graded permutation operator, and the subscript denotes the 4-dimensional auxihary superspace 
V = C2'2 with the grading P\i] = if i = 1, 2 and 1 if i = 3, 4, and 



V = stroK'+{0), 



W = str. 





K4 



(0)i?ft 



(20) 



This impUes that this model, as with the following two model we will study, admits an infinite number of mutually 
commuting conserved currents, thus assuring its integrability. 

The second choice of integrable couplings results from use of an i?-matrix obtained by imposing Z2 grading associated 
with two bosonic and two fcrmionic states to the fundamental su{A) i?-matrix which reads 



/ 



R{u) 



~ 2 
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(21) 



V 

We now solve ( |To| ) and ( pi] ) for K^{u) and Kj^{u). For (pT]), we find that the matrix K^{u) given by 



/ 1 

10 

A^{u) B^{u 

V C_(m) dS{u) j 



(22) 



where 



C-{u) 
D_{u) 



-2u- 4cg + LSajSg + 1) + 1 - 4MSg 

{U + 2Ca - 2.Sa - 1)(U + 2Ca + 2Sa + 1) ' 

^ 

{U + 2Ca - 2Sa - 1){U + 2Ca + 2Sa + 1) ' 

4uS+ 

{U + 2Ca - 2Sa - 1){U + 2Ca + 2Sa + 1) ' 
7/2 - 2w - 4c2 + 4Sa(Sa + 1) + 1 + 4wS^ 



(23) 



[U + 2Ca - 2Sa - 1)(U + 2Ca + 2Sa + 1) 

satisfies (|l0|). The matrix K^{u) can again be obtained from the isomorphism of the superalgebras 71 and 7^ through 

Tf{u)^T_{-u + A). (24) 



Therefore, one choose the boundary matrix K^{u) as 



K+iu) 



(\ 





A+(m) 
yo C+(m) 






D+iu) J 



(25) 
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with 



A+{u) = 



B+{u) 



C+{u) 



D+{u) 



+ 4sa(s6 + l) + 5-4(u-4)Sg 



{u + 2cb - 2sb - 3)(u + 2c6 + 2sb - 1) 

^{u - 4)S,- 

{u + 2cb - 2sb - 3)(u + 2cb + 2sb - 1) ' 



4(w-4)S+ 



{u + 2cb - 2sb - 3)(u + 2cb + 2sb 
P -6u- 4cl - 8cb + 4sa{sb + 1) + 



1)^ 



4(ii-4)Sg 



1) 



(26) 



{u + 2cb — 2sb — 3)(w + 2cb + 2sb 

For this example it can be shown that the Hamiltonian is related to the logarithmic derivative of the transfer 
matrix t(u) with respect to the spectral parameter u at u = (up to an additive chemical potential term) 



L-l 



H^}_ h,.j+i + \k'- (0) 



stroK+{0)HLo 
stroK+{0) ' 



(27) 



with 



2du ^ ' 



and subject to the constraints (^). 

The third choice of integrable couplings results from use of the i?-matrix obtained by imposing Z2 grading to the 
fundamental (7^(3|1) i?-matrix which reads 
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R{u) = 

























u 

















-2 
















-2 

















u 













































-2 








u 






















































-u + 2 






















































—u 








2 


















-2 


























u 




































-2 

















u 













































2 








—u 







v 
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in we solve (|l0|) and 


& 


for K-{u) and if+(u) 


For (§ 


1 we 


obtain 













(28) 



/ 1 \ 

10 

A^{u) B^{u) 

\ C-{u) D_{u) ) 



(29) 



where 



A.{u) 



B^(u) 



C-{u) 



2u ~ Acl + Asa{sa + !) + !- 4uS^ 



{U + 2Ca 



2Sa-l){u- 



2Ca + 2Sa + 1) ' 



{u + 2Ca - 2Sa 



-l){u- 
4uS+ 



2Ca + 2Sa + 1) ' 



D-{u) = - 



{U + 2Ca - 2Sa - 1){U + 2Ca + 2Sa + 1) ' 

-2u- 4cl + ASaiSa + !) + ! + 4mS^ 



(U + 2Ca - 2Sa - 1)(U + 2Ca + 2Sa + 1) 



(30) 
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and 



Tf{u) = jr_(-u, + 2), 



J = diag{\, -1, 1, 1), 



(31) 



giving 



with 



K+{u) 



/ 1 \ 

0-10 

A+{u) B+{u) 

\0 C+iu) D+{u) ) 



(32) 



A+{u) 
B+{u) 
C+{u) 
D+{u) 



-2u- ^cl + Asa{sb + 1) + 1 - 4(m - 2)Sg 
[u + 2ch - 2sf, - 3)(u + 2cb + 2s6 - 1) 

^{u - 2)S,- 

[u + 2cb - 2sb -3){u + 2cb + 2sb - 1) ' 

4(^^-2)8+ 

{u + 2cb - 2sb - 3){u + 2cb + 2sb - 1) ' 
u^-2u- Acl + Asa{sb + 1) + 1 + 4(m - 2)Sg 
{u + 2cb ~ 2sb ~ 3)(u + 2cb + 2sh - 1) 



(33) 



The Hamiltonian (^ is related to the logarithmic derivative of the transfer matrix t(u) with respect to the spectral 
parameter m at w = (up to an additive chemical potential term) 



L-l 1 



5^roX+(G)g^o 
5iroX+(0) ' 



(34) 



with 



For this case we obtain (||) subject to the constraints (||). 



III. THE BETHE ANSATZ SOLUTIONS 



Having established the quantum integrability of the models, let us now diagonalize the Hamiltonians by means of 
the algebraic Bethe ansatz method ]l^ , p7t . For the first case (|^), introduce the 'doubled' monodromy matrix lJ(u) 



U{u) = T{u)K^{u)fiu) 



( A{u) Bi{u) B2{u) Bsiu) \ 

Ci{u) Vii{u) Vi2{u) Vuiu) 

C2iu) V2l{u) V22{u) V23{u) 

\C:i{u) Vsiiu) V32iu) Vsaiu) J 



(35) 



where T{u) = T ^(— u). Substituting into the reflection equation ( [To| ) we may draw the following commutation 
relations, 



Vbd{ui)Bc{u2) 



[Ui - U2 - 2){ui + M2 - 4) 
{Ui - U2){ui +U2-2) 

2(ui - 2)u2 



(Ul+U2-2)(U1-1)(M2-1) 



r(ui +U2- 2)llr{ui - U2)fd^3e{u2)i>g^{Ul) 

ri2ui - 2)f^Bg{ui)A{u2) + 



r{2u,-2)fMu,)Ac{u2), 



A{ui)Bf,{u2) 



iui - U2+ 2){ui + U2) 



-Bp{u2)A{ui) - 



2{ui + U2) 



{Ui - U2){ui + U2 - 2) 



(lil - U2){ui +U2-2) 

2 - 1 

H ■ ■z[Ba{ui){T)ap{u2) -<5a/3^(u2 )] . 

Ul + U2 - 2 ^2 - 1 



Bp{uM{u2) 



(36) 



(37) 
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Here T>bd{u) = t>\jd{u) — :jTrr^6d^(it) and the matrix r{u)^ which in turn satisfies the quantum Yang-Baxter equation, 
takes the form, 



r[[[u) = l, rii{u) = rti{u) 



u + 2 

\i - 2 ' 



u- 2' 



It - 2' 



u-2 



Next choose Bethe state of the form 



(38) 



(39) 



with |0) being the pseudovacuum. Acting the transfer matrix t{u) on the state |f2) we have t{u)\Q) — A{u)\fl) with 
the eigenvalue 



u (cb-f) {cb~ '^ + 2sb + l) -SL {u + Uj){u-Uj+2) 



A{u) 



"-l(cfe-|-l) (c6-f + 2sb) -IjL^ + -2) 



n 



(m — Uj){u + Uj — 2) 



(40) 



provided the parameters {uj} satisfy 



AT 



(cfc - ^) (cb - ^ + 2sb + 1) ^Uj_2^2L ^ -Q (%- - - 2) (mj- + - 4) ^^(1)^^^. i^^i^ 



2(cf,-^-l) (cfc-^ + 2sb) 



(41) 



Here A'^"'^^(w; {u^}) is the eigenvalue of the transfer matrix r(^)(u) for the reduced problem which arises out of the r 
matrices from the first term in the right hand side of ( ^ ) with the reduced boundary K matrices (u) 



(l 







A^_>{u) B'_>{u) 1 , 
\0 C'^\u) d'-^\u) 



(42) 



where 



~ Acl ~ SsgCg + Asq + Aju - l)Sg 

{U - 2Ca){u - 2Ca - ASa ~ 2) 

4(^ - 1)S- 

(u-2Ca)(u-2Ca-4Sa-2)' 

4(»-l)S+ 

{u-2Ca){u-2Ca-4Sa-2)' 

- Acl - 8SaCa + 4Sa - 4(m - l)Sg 
(u-2Ca)(M-2Ca-4s„-2) 



and 



1 

C|'^(u) D^_^\u) 



(43) 



(44) 



where 
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^(1). , ^ -2u- 4cl - 4cb(2sb - 1) + 8sb + 4MSg 
+ ~ {u - 2cb + 2){u - 2cb - 4:Sb) 



{u ~ 2cb + 2){u- 2cb - 4s6) ' 

+ ^ ' {u~2cb + 2){u-2cb-4sb)' 

_ u^^2u~ 44 - 4cbi2sb - 1) + 8sb - 4t.Sg 
+ ^ (u-2cb + 2)(u-2cb-4sfc) ■ ^ ' 

Here /'^^^■'(m), the boundary K matrix after the first nesting, fohows from the relation 
Indeed, applying the monodromy matrix T{u) and its "adjoint" T{u) to the pseudovacuum, we have 



v,bium ^ -^K^Viu)m = K^iuUi—^r'^i^). (46) 



TnHI*) = 1*), TM\^) = ( -)^|*), 

u ~ 2 



fu{u)\^)^0, fdt(u)|*) = 0, fdi(u)|*) = 0. (47) 
where d b, d,b — 2,3, 4. Then we have 

A{um = \^), 

u — 2 u— i u— i 

(u - 1)721 (m)Ti2(w) - T22{u)T22{u) - T^S (w)T32 (u) - r24(u)T42(M) 

= -Tii(u)Tii(m) + (u - l)Tl2(u)T2i(u) - Ti3(u)T3i(w) - fi4(w)r4i(u), 
(li - l)r2i(M)Ti3(w) - T22{u)f23{u) - T23{u)f33{u) ~ T2i{u)fis,(u) = uf'13 (u)T2i (w) 
(li - l)r2i(M)Ti4(w) - T22{u)f2i{u) - T23,{u)fs,i{u) - T2i{u)fii(u) = ufii{u)T2l{u) 
Tzi{u)fi2{u) -{u- l)T32(u)f22(u) + T33(w)f32(u) + T34(u)f42(M) = -MT22 (u)T32 (u) 
Tzi{u)fi3{u) - (-U - l)T32(u)f23(u) + T33{u)f3-i{u) + T34(u)f43(u) 

= T2i(m)Ti2(u) + f22{u)T22{u) + {u + 1)T23 (M)r32 (u) + ^24 ('«)'742 (w) , 
T3i(m)T'i4(u) - (m - l)r32(M)T24(u) + ^33 (m)T34 (u) + r34 (■u)f'44 (w) = uf'24(u)T32 (w) 
T4l(M)ri2(M) + Ta2{u)T22{u) + {u + 1)743 (w)T32 (u) + T44(u)T42(m) = UT32 (u)T43 (m) 
r4l(u)ri3(w) + 742(M)r23(w) + {u + 1)143 (w)T33(u) + 744(^)143(1*) = ufz3{u)Ti3{u) , 
Tii{u)fii{u) + Ti2{u)f2i{u) + {u + l)r43(u)T34(u) + T44(u)T44(m) 

= f3i(M)ri3(u) + f32(u)T23(u) + f33(u)T33(M) + [u + I)f3^{u)n3{u) ■ (49) 

which come from a variant of the (graded) Yang-Baxter algebra (^ with the R matrix (||), 

2 2 

T {u)R{2u) f (u) =f {u)R{2u) T {u). (50) 

Noticing the change u ^ u — 1 with respect to the original problem, one may check that these boundary K matrices 
satisfy the reflection equations for the reduced problem. After some algebra the reduced transfer matrix r'^^ [u) may 
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be recognized as that for the inhomogeneous supersymmetric t — J open chain interacting with the Kondo impurities 
of arbitrary spins, which has been diagonahzed in Ref. [[L3[. The final result is 



AW(u;K}) = 



{cb - f ) (cfc - f + 2sfc + 1)!^ {u-v^ + 2){u + - 2) 



2(cfc-f + 2s6) (c-f-l) AJ-^ + -4) 



n 



M- 1 



N 

n 



Afi 



(u — U.j){u + — 2) T-r (m - + 2)(m + Uq, - 2) 



n 



u — 2 J^-fj- {u — Uj — 2){u + Uj — 4) ^j*-^ (tt — ) (tt + Wq, — 4) 
provided the parameters {vm} satisfy 

Vg (cb - ^)(C{, - ^ + 2sb + 1) (Vg - Uj - 2){Va + Uj - 4) 

Va - I {Cb - ^ + 2sb){cb - ^ - 1) {Va - Uj)(Va + Uj - 2) 



A(2)(^;K},K}) (51) 



= -A(2)(«„;{u,},{«^}). 



(52) 



Here A*^^^(u; {wj}, {wq}) is the eigenvalue of the transfer matrix r*^^'(u) for the M2-site inhomogeneous XXX open 
chain interacting with the Kondo impurities of arbitrary spins, 



A(2)(u;{«^.},{z;4) 



(cb - I) (cb - I + 2sb + 1) 



(cf, - I + 2sf,) (cb 



n 



-7 "T 2 



^ + 2s. 



1) -^-^ c^-f + 2.s^ + l 



Ah 



s " TT {u — wp ~ 'i)(u + wp — i) u — 2 T-r {c^ + ^ — 1) (c^ — -I + 2s^) 
^^7^ j-J-^ (m-w;;3-1)(u + W;3-1) ^ ^j-ft (^7 - I) (c^ + I + - 1) 

TT- (U - I^a)(M + jig - 4) 1^ (m - + 1)(m + W/3 + 1) ^ 

^ iJi + 2)(w + - 2) JJ-^ (m - u;/5 - 1)(u + u;;3 - 1) 

provided the parameters {wp} satisfy 



(53) 



TT (C^ + ^ - ^)(c-y - ^ + 2s-y - ^) l^j- {wp -Va + l)(w/3 + tig ^ 3) _ (wg - - 2){wp +WS -2) 

yilb (^7 - ^ - 5)(c7 + ^ + 2s^ - i) j-Ji («;^ - + 3)(«;^ + - 1) " j-^ («;^ - u;.- + 2){wp + ws + 2) ' ^ ^ 



After a shift of the parameters Uj Uj + 1, + 2, the Bethe ansatz equations (41), ( |52| ) and (p4| ) may be 

rewritten as follows 



n 



1 ^2L "TT (^j - M» + 2)(U.,- + + 2) _ {Uj - fg + l)(Mj + fg + 1) 

1 ^J- {uj -Ui- 2){uj +Ui -2) {uj - Va - 1)(% +Va-l)' 

+ ^ + 2s^ 1^ (Ug - + 1)(Vq + + 1) _ y-f (Wg - + l)(Wg + + 1) 



I. ^7 



n 



n 



^ + 2S^ -^^-^ (Ug - - l)(Wg + Uj - 1) ^J-^ (Wg " W/5 " l)(Ug + tW/J " 1 ) ' 

n (c^ + ^ - \) (c^ - ^ + 2s^ - ^) ^ {wp -Vg-l) {wp + i;g - 1) _ TT" [wp-ws- 2) [wp + ws-2) 



) ^J-^ (?«/3 - I'g + 1) (wp + -Wg + 1) ^^-^ (w/3 - W5 + 2) {wp + W5 + 2) ■ 



(55) 



with the corresponding energy eigenvalue E of the model 

N 



E 



1 "'-I' 



(56) 



We now perform the algebraic Bethe ansatz method [^2|jT^] procedure for the second couplings (||). Introducing 
the 'doubled' monodromy matrix U{u), 



U{u) = T{u)K^{u)f{u) = 



( A{u) Bi{u) B2{u) B^iu) \ 

Ci{u) Viiiu) Vuiu) Vu{u) 

C2{u) V2l{u) V22{u) V23{u) 

\C3iu) Vsliu) V32{u) 2?33(u) / 



(57) 
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where T{u) = T ^(—w). Substituting into the reflection equation (p^,we may draw the following commutation rela- 
tions, 



(U1+U2-2)(U1-1)(U2-1) 

-r{2ui-2)f^Bg{ui)V,,{u2), (58) 



(m1 - U2)(mi - 1) 

A( \K I \ - U2 + 2){ui + U2) , , , , 2(ui+U2) 

■^(W1)S/3(W2) = 7 77 —B0{u2)A{ui) - —Bp{ui)A{u2) 

[Ui — U2)(Ui +U2 — 2) (Ul — U2)(Ui +U2 - 2) 

+ — ::[Bc,{ui)if>af3iu2) ^Sc,pA{u2)]. (59) 

Ul + U2 - 2 U2 ~ I 

Here I?bd(w) = Pfcji(u) — :;;^SbdA{u) and the matrix r{u) ,which in turn satisfies the quantum Yang-Baxter equation, 
takes the form, 



rlliu) 






1 — ' 








= r\l{u) 


= rlliu) 


= rlliu) 


= rf,iu) 


- ^32 

— '32 




= rl\{u) 


^ rlliu) 


^ rlliu) 


-rlliu) 


— '23 



u-2' 
u 



(60) 



u-2 

Choosing the Bethe state as 

\n)=BM---B,„iuN)\Q)F'^-'«, (61) 

with |0) being the pseudovacuum, and acting the transfer matrix t(u) on the state |r2),we have t{u)\VL) = A(u)|ri),with 
the eigenvalue, 

K( \ - " - 4 (cfc + f + Sb + |) (cb + f - Sb - ^) 1^ (m + Uj)iu - Uj + 2) 
^ (cb + f + Sb-i) (cb + t - ,Sb - I) (u-7.j)(" + %-2) 

u — 1 u — 2 -"^^ (u — Uj)(u + — 2) 

provided the parameters {wj} satisfy 

- 4 (cb + ^ + sb + I) (cb + ^ - Sb - \) Uj - 2 2L _ -rr ("j ~ ^ 2) (uj + - 4) ^ ^ 
u,-2{cb + ^ + Sb-\)icb + ^-Sb-^)^ u, ' l\iu,-u, + 2) iu,+u,) ' ^^"^^ 

Here A(^'(u; {ui}) is the eigenvalue of the transfer matrix T(^'(7i) for the reduced problem, which arises out of the r 
matrices from the first term in the right hand side of (^),with the reduced boundary K matrices K^^\u) as, 

/ 1 

A^}\u) B^_}\u) I , (64) 



K^l\u) 



\o c):>iu) D'y{u) 



where 



.(1), , u" -Ac^^~ACa+^SaiSa + l)+i-^iu-l)^l 

A_ [u) — — 



iu + 2Ca - 2Sa - 1)(W + 2Ca + 2Sa + 1) 

4(^ - 1)S- 

iu + 2Ca - 2Sa - 1)(U + 2Ca + 2Sa + 1) ' 

4(^^-1)8+ 

iu + 2Ca - 2Sa - l)iu + 2Ca + 2Sa + 1) ' 



^(1). ^ _ »'-4cg-4Ca+4ga(s, + l)+3 + 4(M-l)Sg 
iu + 2Ca - 2Sa ~ l)iu + 2Ca + 2Sa + 1) 
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and 



where 



1 

K^^\u)= I A^_^\u) B^_^\u) I , (66) 
d+\u) D^'^(u) 



6u - Acl - 8cb + 4sfa(s6 + 1) + 5 - 4(m - 4)Sg 



(« + 2cb - 2sb - 3)(u + 2cb + 2sb - 1) 



(u + 2c6 - 2sb - 3)(« + 2c6 + 2sb - 1) 



{u + 2cb - 2sb - 3)(m + 2cb + 2sb - 1) ' 

+ (m + 2c6-2s6-3)(u + 2c6 + 2s6-1) ' ^ ^ 

Here ii'i^^(u), the boundary K matrices after the first nesting, follows from the relations, 

Vd,{um = -^K'^^ium = K.iuUi-^f'l^)- (68) 

ct -L U Zi 

Indeed, applying the monodromy matrix T(u) and its "adjoint" T{u) to the pseudovacuum, we have 

Tn(«)|*) = 1^), T,,ium = i^)^m, 

u — 2 

TUum + 0, rdb(«)|*) = 0, Trfi(u)|*) = 0, 
fn(u)|*) = |vl/), = (-^)^|*), 

Tid(w)|*)^0, Tk6(u)|*)=0, Tdi(u)|*)=0. (69) 
where d^h, d,b = 2,3, 4. Then we have 

^(^)|*) = I*), 

Sd(u)|*)^0, Cd{u)\^)=0, 

■DM^m = (^)^^i^-(«)<i6i*), 

'D.dium = {-^r'^iK_iuU+-^)\^) - -^1*). (70) 

U ^ U J. U 1. 

{U - l)T2i(u)Ti2(u) - T22(M)r22(w) - T23{u)fs2{u) - T24(w)T42(m) 

= -fniu)Tn{u) + {u - l)Ti2(u)T2i(u) - T'i3(u)T3i(u) - Ti4(u)T4i(u), 

{U - l)T2l{u)Ti3{u) - T22{u)T2z{u) - T23{u)T33{u) - T2A{u)Ti'i{u) = uTi3{u)T2i{u) 
{U - l)T2l{u)fi4{u) - T22(u)f2i{u) - T^iS (w)?34 («) - T24 (m)T44 (w) = uTm (m)T21 (u) 
T3i{u)fi2{u) - {U- l)T32{u)f22{u) + (w)T32 (u) + T34(m)T42(m) = -WT22 (m)T32 (u) 
T3l(w)Tl3(«) - (« - l)T32(u)f23(u) + T33(«)T33(w) + T34(w)f43(«) 

= f2i(u)ri2(u) + T22(u)r22(«) - {u - l)f'23(u)T32(«) + T24(m)T42(u), 
T3i(u)Ti4(u) -{U- l)T32{u)f2i{u) + T33{u)f3i{u) + T3i{u)fii{u) = -uf2i{u)T32{u) 
Tii{u)fi2{u) + T42(m)T'22(m) - (u - l)r43(w)f32(M) + T44(u)'f42(w) = ^7132 (u)T43 (u) 
Tii{u)fi3{u) + T42(M)f'23(u) - (m - l)T'43(M)T'33(?i) + T44 (71)743 (it) = -"Tbs (w)l43 (") , 
T4i(M)ri4(u) + r42(M)T'24(M) - (" - 1)T'43 (w)^34 («) + T44(?i)T44(u) 

= f3i(M)Ti3(u) + T32(m)T23(u) + T33(u)r33(M) - {u - l)f3A{u)Ti3{u) . (71) 
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which come from a variant of the (graded) Yang-Baxter algebra with the R matrix (^l|), 

2 2 

T {u)R{2u) f (u) =f {u)R{2u) T (u). (72) 

Noticing the change u ^ u — 1 with respect to the original problem, one may check that these boundary K matrices 
satisiy the reflection equations for the reduced problem. After some algebra, the reduced transfer matrix t^^^{u) may 
be recognized as that for the inhomogeneous su(3) t — J open chain interacting with the Kondo impurities of arbitrary 
spins, which has been diagonalized in Ref. The final result is, 

^(1), . r |^ ^ U - 4 (Cfc + f + Sb + ^) (Cb + I - Sb - ^) ^ (u - Vc, + 2){u + - 2) 
[^A^ji) y_2(cb + f +Sb-i)(cb + |-Sb-|) J-j-^ {U - Va,)iu + Vc. - A) 

U — 2 (u — Uj —2)[U + Uj — 4j (u — Wa) (u + ■t'a — 4) 

provided the parameters {wm} satisfy 

- 4 (Cb + ^ + 3fc + ^) (Cfc + ^ - Sb - \) 1^ (Vg - Uj - 2){Va + Uj - 4) -p- {Vg - VC, + '^){Va + ~ ^) 

-3 (cb + ^ + Sb - \) {cb + ^-Si- I) K +Mj -2) I i K-wc-2)(w„+v^-6) 

= -A(2)(z;„;{«J,{«4). (74) 

Here A*^^)(u; {wj}, {wq}) is the eigenvalue of the transfer matrix t^'^^{u) for the M2-site inhomogeneous XXX open 
chain interacting with the Kondo impurities of arbitrary spins, 

A(2)r,,-/„ \ {„ 1 ^ - i^b + ^ + Sb + ^) {Cb + ^ - Sb - ^) yr C^, - § + + | 



M. 



_ 1 



u - A YT {u - Wfj + 2){u + - A) u-2 -p-r (c^ + § + - |) (c-y - | - s-y + |) 

^M-3 jj-^ {U-Wp){u + Wp-Q) ^ yiab (C7 + f - S7 - 5) (Cy - | + Sy + |) 

^ TT (u-«a)(M + fa -4) |^ (m ~ - 2) (u + - 8) ^^^^ 

(u - Ua - 2)(m + - 6) J-J-^ (u - W/3)(u + Ui/S - 6) 
provided the parameters {wp} satisfy 

n (Cy + ^ + - |)(Cy - ^ - S-y + I) jT- (Wg -i^a)(w/3 +t;a -4) ^ -TT ("'/3 ^ ™<5 + 2)(W/5 + - 4) 



via b (C7 - ^ + + f )(c7 + - •'^7 - 5) (^/S - - 2)(w^ + - 6) (u;/3 -ws- 2){wfs + 



(76) 



After a shift of the parameters Uj ^ + 1, Vm — > Wm + 2, ^ + 3, the Bethe ansatz equations ( |63| ) , ( |74| ) and 
( [76| ) may be rewritten as follows 

^ ^J ~ 1 ^2L TT (Mj " Mi + 2)(Uj + Uj + 2) _ -pj- {Uj -Va + l)(Uj +^7^ + 1) 
+ 1 (Wj - - 2){u.j + Ui - 2) (Uj - - l)(uj + - 1) ' 

y-r C^, + ^ + S^, + I -jX- (t;^ - - l)(fa + Uj - 1) _ |t- (fg - W/3 + l)(z7a + W/j + 1) 
4i"b C7 - ¥ + «7 + 1 y K - + + Uj + 1) ~ (Wa " W/3 " l)(Ua + - 1) ' 



^ TT (^g - - 2)(t;a +VC-2) 
(ua -«;<: + 2)(wa + + 2) 



n (cT, + ^ + + l)(c^ - ^ - + 1) 1^ (w/3 - Ug + 1) {wp + Vg + 1) _ 1^ {wp -ws + 2) {wf3 +WS + 2) 

-y=a,b - ^ + ''^7 + 1)(C7 + ^ - Sy + 1) (u;^ " l^g " 1) (u^^ + fg - 1) " \}^ {w p - Ws ~ 2) (w0 + U-5 - 2) ' ^ ^ 
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with the corresponding energy eigenvalue E of the model 



N 



^ M? - 1 



(78) 



1 J 



We now perform the algebraic Bethe ansatz method procedure for the third couplings (^). Introducing the 

'doubled' monodromy matrix U{u), 



U{u) = T{u)K^{u)T{u) 



( A{u) Bi{u) B2{u) B^iu) \ 
Ci{u) Vn{u) Vu{u) Vuiu) 

C2{U) V2l{u) V22{U) I?23(") 



(79) 



where T(u) — T ^(—m). Substituting into the reflection equation (p^,we may draw the following commutation rela- 
tions, 



iibd{ui)Bc(u2) = 



(Ul - U2 - 2)(mi + U2) 
{Ul — U2){ui + U2 + 2) 
2uiU2 



(Ul +U2 +2)(ui + 1)(m2 + 1) 



riu, +U2 + 2)lir{u, - U2)fdBe{u2)VAui) 
r(2^il+2)fXKM(^2) + 



A{ui)Bp{u2) 



{Ui - U2){ui + 1) 

(ui -~U2- 2)(wi + U2) 



r{2ui + 2)llBg{ui)V,,{u2), 



Bp{u2)A{ui] 



2{ui + U2) 



[ui - U2){ui + U2 + 2) 

2 

\ (mi ) (2?a/3 (W2 ) 

Ml + M2 + 2 



(ui - U2){ui + U2 + 2) 



Bi3{ui)A{u2) 



U2 — 1 



5al}A{u2)\. 



(80) 



(81) 



Here I'm(u) = Pft(i(u) + ■;^^5i,4A{u) and the matrix r{u) , which in turn satisfies the quantum Yang-Baxter equation, 
takes the form, 



r\\{u) 


= rll{u) 


= rll{u) 


= 1, 










^rl\{u) 






- ^32 

— '32 


r\l{u) 


= rl\{u) 


= rll{u) 




= rll{u) 


_ ™,32 
^ '23 



u-2' 
u 



u-2 



Choosing the Bethe state |0) as 



\n)^B,,{ui)---B,^{uNmF''- 



(82) 



(83) 



with |0) being the pseudovacuum, and acting the transfer matrix t[u) on the state |r2),we have t{u)\^) = A(u)|f2),with 
the eigenvalue, 



2{cb + ^-Sb-\){cb + ^ + Sb + \) ^{u + Uj)(u - uj - 2) 



u+l{cb + ^ + Sb-^) (cb + ^-Sb-^) {u - Uj){u + Uj + 2) 



N 



\2L 



J|(!i±M!i_^^A«(M;K}), 



u + 2 [u — Uj){u -\- Uj -f- 2) 



(84) 



provided the parameters {uj} satisfy 

Uj - 2 (cb + ^ - Sfc - i) (cfc + ^ + Sfc + i) Uj + 2 



11 

u, {cb+'^ + Sb-\){cb + ^-Sb-^' 



i- 



Here A'^^(u; {ui}) is the eigenvalue of the transfer matrix r'^'(u) for the reduced problem, which arises out of the r 
matrices from the first term in the right hand side of (pO|),with the reduced boundary K matrices k2:\u) as. 
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1 

(1)^ 



K^^iu) = I A'_!'{u) Bl''{u) I , 







d^iu) D'y{u) 



where 



and 



where 



- 4cl + 4Ca + 45a (g, + 1) ^ 1 ~ A{u + 1)S^ 

{u + 2c„ - 2.s„ - 1) (?i + 2c„ + 2sa + 1) 
^(" + l)S,7 

(U + 2Ca - 2Sa - 1){U + 2Ca + 2s„ + 1) ' 

4(^^+1)8+ 

(it + 2Ca - 2Sa - 1)(W + 2Ca + 2Sa + 1) ' 

- 4cl + 4ca + 4s,{sa + 1) - 1 + 4(u + l)Sg 

{u + 2Ca-2Sa-l){u + 2Ca + 2Sa + l) 



/-I 



\ 

(1)^ 



av^(m) 

V ci^^H D^^\u)J 



^(1). . ^ - 2m - 4cg + Asbjsb + 1) + 1 - 4(m - 2)Sg 
+ {u + 2cb - 2sb - 3)(u + 2c6 + 2sb - 1) ^ 



A{u ~ 2)S- 



{u + 2cb - 2sb - 3)(u + 2cb + 2sb - 1) ' 

4(t^-2)S+ 

{u + 2cb - 2sb - 3)(u + 2cb + 2sb - 1) ' 



^(1). . ^ u''-2u- 4cl + 4.%{sb + l) + l + 4{u- 2)Sg 
+ ^"'^ {u + 2cb - 2sb - i){u + 2cb + 2sb - 1) 

Here ii'i^^u), the boundary K matrices after the first nesting, follows from the relations, 



^db(w)|*) 



u+1 



dd 



-)( 



u+l u+2 



Indeed, applying the monodromy matrix T{u) and its "adjoint" T{u) to the pseudovacuum, 

Tid(w,)|*) ^ 0, rrffc(u)|*) = 0, Tdi{u)\^) = 0, 



fnHI*) = I*), fad{um - (- 



u + 2' 



TidMI*) 7^ 0, Td(,(tt)|*) =0, Tdi(M)|*)=0. 
where b, d,b = 2,3, 4. Then we have 



Vab{u)\^) = (- 



u + 2 



Vdd{um = {-^^r{K_iu)dd-^^^ 



u+l 



If). 
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{u + l)T2l{u)fi2{u) + T22{u)f22{u) + T2z{u)f^2{u 

= fii{u)Tii{u) - (m - l)Ti2(u)T2i 

{u + l)r2i(li)Ti3(M) + T22{u)f23{u) + (-u)f'33 (u 
{u + l)T2l{u)fi4(u) + T22{u)f2i{u) + T23iu)f34{u 
Tzi{u)fi2{u) - (m - l)T32(M)f22(u) + T'33(m)T32(u 
?3l(w)Tl3(u) - (m - 1)732 (u)T23(u) + r33(w)T33(u 



T2,i{u)Ti4{u) - {u - l)T32(u)T24(u) + T33{u)T34{u 
T4i(m)Ti2(m) + T42(u)f'22(M) - {u ~ l)r43 (m)T32 (u 
T4i(m)Ti3(u) + T42(u)T'23(w) - (u - 1)T43 (m)T33 (u 
T4i(u)Ti4(w) + T42(u)f'24(w) - - 1)T'43 (u)T34 (u 

= T3i(m)Ti3(u) + f32{u)T23[u) + T33 (u)T'33 (") - (u - l)f'34(u)T43(u). 

which come from a variant of the (graded) Yang-Baxter algebra with the R matrix (p8|), 



+ T24(u)T42(w) 
u) + f'i3(M)r3i(u) + f'i4(w)r4i(u), 
+ r24(u)T'43(M) = -ufi3{u)T2l{u) 

+ T2i{u)fii{u) = -ufu{u)T2i{u) 

+ T3i{u)fi2{u) = -uf22{u)T32{u) 
+ T34(u)f43(M) 



T2i(u)Ti2(u) + r22(u)T22(w) - {u ~ 1)T23 («)r32 (u) + T24(u)r42H, 



+ T34(u)T44(m) = -mT24(m)T32(u) 
+ T44(u)T42(m) = uf32{u)Ti3{u) 
+ T44(u)T43(u) = -mT33(m)T43(u), 
+ T44(u)r44(u) 



(93) 



T {u)R{2u) f [u) =f {u)R{2u) T (u). 



(94) 



Noticing the change u — » u + 1 with respect to the original problem, one may check that these boundary K matrices 
satisfy the reflection equations for the reduced problem. After some algebra, the reduced transfer matrix t^^^{u) may 
be recognized as that for the inhomogeneous su(3) t — J open chain interacting with the Kondo impurities of arbitrary 
spins, which has been diagonalized in Ref. [O]. The final result is. 



aW(«;K})- 



2 (cfe + f - Sfe - i) (ch + f + Sb + 5) -i^ (u - w„ + 2)(m + Va + 2) 



(C5 



Sb 



n 



2) a=l 



U 



U 



1 ^ 

-n 



^ (u — Uj){u + Uj + 2) (u 



{u ~ Va){u - 

2){u + Va - 2) 



J- J; {U-Uj - 2){U + Uj) ^\ 

provided the parameters satisfy 

N 



A(2)(u;K},K}) 



(95) 



Vn 



2 (cb + - Sb - i) (cb + ^ 



Sb 



Vn 



l(Cb 



Sb 



n 



2){Va + Uj) {Va " t^f + 2){Va + + 2) 

Vc;-2){Va + V^-2) 



n 



= -A^'\vn,;{u.,},{vp}). 



(96) 



Here A^'^\u; {uj}, {va}) is the eigenvalue of the transfer matrix r*^^-'(u) for the M2-site inhomogeneous XXX open 
chain interacting with the Kondo impurities of arbitrary spins, 
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After a shift of the parameters Uj 
rewritten as follows 



l,w„ 



Wm + 1, the Bethe ansatz equations (|96|) and (Bq) may be 
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with the corresponding energy eigenvalue E of the model 
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IV. CONCLUSION 

In conclusion, we have studied integrable Kondo problems describing two boundary impurities coupled to one- 
dimensional extended Hubbard open chains. The quantum integrability of these systems follows from the fact that 
the Hamiltonians in each case are derived from a one-parameter family of commuting transfer matrices. Moreover, 
the Bethe Ansatz equations and expressions for the energies are derived by means of the algebraic Bethe ansatz 
approach. We would like to emphasize that the boundary K matices found here are non-regular in that they can not 
be factorized into the product of a c-number K matrix and the local momodromy matrices. However, similar to the 
cases discussed in [^,^, it is possible to introduce a singular local monodromy matrix L{u) to express the boundary 
K matrix (u) as. 



K^{u) = L(u)Z~i(-u). 
where, for example in the case of the superalgebras gl{2\2) model 
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(102) 



1 - 2S^ / 



which constitutes a realization of the Yang-Baxter algebra (^ when e tends to 0. The recent work of Frahm and 
Slavnov JT^ confirms the existence of such non-regular solutions by means of a projecting method. 

Finally, we would like to stress that here we have only considered the case of Kondo impurities in these extended 
Hubbard models which are based on the sl{2) subalgebra of the bulk symmetry of the models. It is of course 
possible to consider other boundary impurities corresponding to different subalgebra embeddings such as sZ(l|l) for 
the gl{2\2), gl{3\l) cases or sl{3) for the gl{3\l), gl{4:) models and even gl{2\l) for gl{3\l), gl{2\2). For the case of 
t ~ J models such other types of integrable boundary impurities have been studied in . 
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APPENDIX A: DERIVATION OF THE NON-C-NUMBER BOUNDARY K-MATRICES 



In this appendix, we sketch the procedure of solving the (Z2-graded) RE for K^{u) . To describe integrable Kondo 
impurites coupled with the one-dimensional supersymmetric extended Hubbard model open chain, it is reasonable to 
assume that 



/ 1 \ 

10 

A{u) B{u) 

V C(m) D{u) J 



(Al) 



Throughout, we have omitted all the subscrips for brevity, reflecting that the fermionic degrees of freedom do not 
occur, as it should be for a magnetic impurity. For the i?- matrix (||), one may get from the RE ( p^ ) 54 functional 
equations, of which 14 are identities. After some algebraic analysis, together with the su{2) symmetry, we may assume 
that 



+ /3(u)S^ B{u) = P{u)^- , 
C(w) = /3(m)S+, =a(u) -/3(m)S". 



(A2) 



There are 10 equations automatically satisfied and 10 same equations , leaving only 20 equations left to be solved 
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- u+u.{C{ui)D{u2) - D{u2)C{ui)), 



2u^{A{u2)A{ui) 
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- D{ui)C{u2)) = 2u+{A{u2)C{ui) - A{ui)C{u2)) - 

- D{u2)C{ui)) = 2u+{C{ui)D{u2) - C{u2)D{ui)) 

- B{u2)C{ui) - C{ui)B{u2) - D{ui)D{u2)) 

2u+{A{ui)D{u2) - A{u2)D{ui)) - u+u^{B{u2)C{ui) - C{ui)B{u2)), 

- B{ui)C{u2) - C{u2)B{ui) - D{u2)D{ui)) 

= 2u+{D{u2)A{ui) - D{ui)A{u2)) - u+u_{B{ui)C{u2) - C{u2)B{ui)), 
2u-{A{ui)B{u2) + B{ui)D{u2)) + u+u^{A{ui)B{u2) ~ B{u2)A{ui)) 
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2u^{A{u2)B{ui) + B{u2)D{ui)) + u+u^{B{u2)D{ui) - D{ui)B{u2)) 
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- 2u+{Diu2)C{ui) - D{ui)C{u2)) + A{C{u2)A{ui) + D{u2)C{ui) - C{ui)A{u2) ~ D{ui)C{u2)), 
{C{u2)A{ui) + D{u2)C{ui)) + u+u^{C{u2)A{ui) - A{ui)C{u2)) 

- 2u+{C{ui)A{u2) - C{u2)A{ui)) + A{C{ui)A{u2) + D{ui)C{u2) - C{u2)A{ui) - D{u2)C{ui)), (A3) 



2u( 



2u 



with w+ = ui + U2, u_ = Ml — U2. Substituting (A2) into these equations, we find that all these equations are reduced 
to the following three equations 
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u+(a{ui) - a{u2)) = + a{ui)a{u2) + s{s + 1)/3(mi)/3(m2)), 

u+{l3{ui) - /3(u2)) = u-(a(ui)/3(M2) + a(u2)/3(Mi) - /3(ui)/3(u2)), 
2u+{a{u2)(3{ui) - a{ui)P{u2)) = 2u_(a(ui)^(M2) + a(M2)/3(ui)) - m_(u+ + 2)/3(ui)/3(w2)). (A4) 



Taking the limit ui 1*2, these equations become 



da{u) 


1 


du 


^ 2u 


d(3{u) 


1 


du 


~ 2u 


da{u) 


1 


du 


~ 251 



{-l + a{uf + s{s + l)l3{uf), 



a{u)^^ - (3{u)^^^ = ^(2a(M)/3(u) - + (A5) 



Solving the first two equations, we have 



, , (ciC2 - M^)(2s + 1) + (C2 - Ci)u 2(C2 - Ci)u 
Q!(w) = 77: 77 77 ^ , P[u) = 77: 77 77 7, (Ad) 

^' (2s + l)(ci -m)(c2-u) ^' (2s + l)(ci -u)(c2-u) ^ ' 



where ci and C2 are integration constants. Substituting these Results into the third equation in (A5),we may establish 
a relation between ci and C2: C2 = ci — 4s — 2. This is nothing but the non-c- number boundary K matrix ( p^ (after 
a redefinition of the constant: ci — > 2c + 4s + 2). 

A similar construction also works for the quantum R matrix and (p8|). 
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